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We analyze the electronic properties of a two-dimensional electron gas rolled-up into a nanotube 
by both numerical and analytical techniques. The nature and the energy dispersion of the electronic 
quantum states strongly depend upon the geometric parameters of the nanotube: the typical radius 
of curvature and the number of windings. The effect of the curvature results in the appearance of 
atomic-like bound states localized near the points of maximum curvature. For a two-dimensional 
sheet rolled up into an Archimedean spiral we find that the number of bound states is equal to the 
number of windings of the spiral. 
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I. INTRODUCTION 

As a consequence of relaxation of elastic stresses, a thin 
solid film that is subject to compressive strain, curls up 
after being partially released from its substrate^. This 
occurs when the strain gradient across the film thickness 
is sufficiently large to overcome the tendency to form 
wrinkles, which appear in the opposite limit of small 
strain gradient 3 . The roll- up of a thin solid film into a 
rolled-up nanotube (RUNT) is particularly exciting since 
RUNTs have a unique structure^^ that mimic the cylin- 
drical symmetry of a radial crystal. This is reflected 
in their peculiar magnetic 6 9 and optical properties. 
Moreover, RUNTs are promising candidates for applica- 
tions in fields ranging from nanofluidics to optics^^. 
The experimental progress in manufacturing rolled-up 
nanostructures triggers the need for a comprehensive the- 
oretical understanding of the quantum carrier dynamics 
in curved nanomaterials. 

The formal description of the quantum motion of a par- 
ticle confined to a curved surface was a puzzle for a long 
time. The problem arises because Dirac quantization in a 
curved manifold leads to operator ordering ambiguities^. 
The situation was cleared up by Da Costa in Ref. 16. The 
formal description becomes well-defined when the con- 
finement of the particle on a curved n-dimensional mani- 
fold is treated as the limiting case of a particle inan + 1 
dimensional manifold that has a confining force acting 
in the normal direction of its n-dimensional surface. Be- 
cause of the lateral confinement, quantum excitation en- 
ergies in the normal direction become much higher than 
in the tangential direction. Henceforth, one can safely 
ignore the particle motion in the direction normal to the 
surface and on the basis of this deduce an effective, di- 
mensionally reduced Schrodinger equation. This proce- 
dure is obviously the most rigorous and physically sound 
one for two-dimensional (2D) curved systems embedded 
in an ordinary euclidean three-dimensional space. In this 
case one finds that due to the curvature a scalar poten- 



tial of purely quantum nature appears in the effective 2D 
Hamiltonian. Its magnitude is related to the local sur- 
face curvature^ so that the quantum mechanics of parti- 
cles confined to thin curved layers is different from those 
on a flat plane. Several studies have analyzed the influ- 
ence of the curvature induced Da Costa scalar potential 
on the electronic states^^S] and the electron transport 
properties'^ of a number of different curved systems 
with complex geometrical shapes. Particularly interest- 
ing is the interplay of curvature and electron-electron in- 
teraction effects 27 . 

Here we concentrate on rolled-up nanostructures, in 
particular in the form of Archimedean spirals. Al- 
though single material structures have been proposed^ 
and even fabricated 2 ^, RUNTs are generally made from 
bilayer or multilayer thin films of different materials, 
e.g. GaAs/InGaAs. The two-dimensional electron gas 
(2DEG) in one the layers is thus confined on a cylindri- 
cal surface whose cross section can be fairly approximated 
by an Archimedean spiral r = l<p where r and <fi are the 
cylindrical coordinates in the plane perpendicular to the 
cylinder axis z and I is related to the radial superlattice 
constant by a r = 2ttI, see Fig.JjJb). The aim of this work 
is to investigate the single particle states of a 2DEG in a 
RUNT. The characteristic Coulomb-like form of the cur- 
vature induced scalar potential implies the appearance 
of localized, atomic-like states. We investigate how their 
corresponding binding energies are related to the length, 
curvature, and inner radius of the nanotube and proof 
that the number of these bound states is equal to the 
number of windings of the spiral. 

This paper is organized as follows: in Sec. [H] we in- 
troduce the geometry of the system under study and the 



correspondent effective Hamiltonian; in Sec. Ill the the- 
ory is applied to calculate spectra and wavefunctions and 
we conclude in Sec. HVl 
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II. HAMILTONIAN OF A 2DEG IN A RUNT 

We first derive the effective Hamiltonian for electrons 
bound to the surface S of a RUNT. As discussed in 
the previous section, electrons in a RUNT are confined 
to a cylindrical surface whose cross section can be ap- 
proximated by an Archimedean spiral [see Fig. [TJb)] . 
Therefore it is natural to adopt cylindrical coordinates 
r = {r, (/>, z} and parametrize the surface S as 
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with z G (—oo, oo ) whereas <\> G (0i n , <\>out)- The endpoint 
of the Archimedean spiral <j)i n (4>out) is related to the in- 
ner (outer) radius of the RUNT by Ri n ,out — l<t>in,out 
where I is the typical length scale of the radial super- 
lattice constant a r = 2ttI . The maximum radius of the 
outer tube rotation is instead related to the number of 
rotations Nr by 



R 



R in + 27rlN R , 



where Nr is treated, for convenience, as a continuous 
variable. From Eq. 0, the covariant components of the 
surface metric tensor are 



g^ = l 2 (l + </> 2 ) 

9z,z 1, 

gcf>,z = g z ,(f> = o, 
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whereas the covariant components of the Weingarten cur- 
vature tensor^ come out 
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FIG. 1. (Color online) (a) Sketch of the surface S of a RUNT. 
The Cartesian reference frame we use throughout this paper 
is indicated, (b) Cross section of the RUNT perpendicular to 
the cylindrical axis. Ri n , R ou t correspond to the inner and 
the outer tube rotation radius respectively. 



where JC is the kinetic energy operator for a particle 
constrained to move along a planar Archimedean spiral 
waveguide 
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whereas Vg is the attractive geometric potential induced 
by the curvature^ 
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One should note that the kinetic energy term Eq. (|6| 
and the geometric potential Eq. ^ are different from the 
expressions derived previously^ by non-trivial numerical 
factors. In the following sections, we will find the eigen- 
states of the Hamiltonian Eq. ([5| by imposing on the i/j(<f>) 
component of the surface wavefunction Dirichlet bound- 
ary conditions at the inner and outer radius of the RUNT 
and requiring, as usual, square integrability. 



The mean curvature is then given by M = a^^/2 
whereas the Gaussian curvature is obviously zero. Fol- 
lowing Ref. [3TJ the effective 2D Hamiltonian for the tan- 
gential motion to the surface S becomes: 
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where m is the effective mass. Since the translational 
invar iance along z remains unbroken, the surface wave- 
function separates as 



ik z z 
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where k z is the momentum along the RUNT axis. This 
leads to an effective one-dimensional (ID) Hamiltonian 
for the component of the surface wavefunction 



2^2 



H 1D = /C + V G 
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III. CURVATURE INDUCED BOUND STATES 

The bandstructure corresponding to the effective ID 
Hamiltonian for a 2DEG in a RUNT in Eq. ([5| consists 
of parabolic subbands 



E n (k z ) = E° n 



2 7,2 



2m 



with n denoting an integer subband index (n > 0). Note 
that the zero of the energy has been fixed at the bottom 
of the 2DEG conduction band in its planar configura- 
tion. It is then obvious that our problem reduces to the 
motion of a particle along a planar Archimedean spiral 
where the subband index n and E® respectively label 
the eigenmodes and the corresponding eigenenergies of 
the Hamiltonian: 



H° = /C + V G . 



(8) 
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FIG. 2. (Color online). (a) Electronic spectrum of an 
Archimedean spiral with inner radius Ri n = 3tvI and differ- 
ent numbers of rotations. The points are the numerical re- 
sults of the exact diagonalization of the Hamiltonian whereas 
the corresponding continuous lines are the energies obtained 
considering the effect of the geometric potential in first-order 
perturbation theory, (b) Same for a RUNT with inner radius 
Rin — 32 ttL 



The exact eigenstates of H can be found by writing the 
ip(<fi) component of the total surface wavefunction as 



i>(<t>) 



(9) 



where the \j ' s are the eigenstates of the kinetic energy 
operator. To proceed further, it is convenient to intro- 
duce the arclength of the Archimedean spiral measured 
from 6 = 



3(0) 



<j>(l + tf) + log (<f> + y/l + </> 2 ) 



(10) 



In terms of s the kinetic energy operator takes the com- 
pact form /C = —h 2 d 2 / (2m) and the corresponding eigen- 
states can be written as standing waves 
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In the equation above L indicates the total length of the 
Archimedean spiral whereas Si n is the arclength value at 
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FIG. 3. (Color online). Comparison of the exact (full line) 
and the asymptotic (dashed line) form of the geometric po- 
tential measured in units of h 2 /(2ml 2 ) as a function of s/l. 
Inset: same for the behavior of the arclength as a function of 
the azimuthal angle <j> of the spiral. 



the inner radius of the RUNT. By direct diagonalization 
of the Hamiltonian on the basis of the Xj' s ? we obtain 
the eigenstates and the corresponding energy spectrum 
for any value of Ri n and Nr, choosing I as the unit length 
scale. All reported calculations are obtained introducing 

a CUtoff jmax 

100 in the infinite sum Eq. ([9|, which in 
all cases is sufficient for convergence. 

As shown in Fig. [2] the spectrum consists of two dis- 
tinct regions. At high energies, the spectrum has a free- 
particle-like quadratic dependence on n (E® ~ n 2 ). In 
this regime, a good approximation consists in retaining 
the effect of Vg in first-order perturbation theory [contin- 
uous lines in Fig. [2]. On the contrary, the low energy part 
of the spectrum is dominated by the effect of the geomet- 
ric potential which therefore produces a strong mixing of 
the free particle states. For integer number of rotations, 
there is a critical mode that separates the two extreme 
spectral structures corresponding to n = Nr indepen- 
dent of the inner radius of the RUNT. This critical state 
corresponds to a zero energy state where the geometric 
potential energy balance the kinetic energy. 

Next we show that the appearance of these two dis- 
tinct spectral structures emerges as a natural conse- 
quence of the competition between the confinement due 
to the Dirichlet boundary conditions at the inner and 
outer radius of the RUNT and the effect of the geomet- 
ric potential. In order to make a qualitative analysis of 
the spectrum, it is convenient to consider the asymptotic 
form of the geometric potential Vg ~ —h 2 /(Sml 2 (j) 2 ). 
Apart from a logarithmic correction, the arclength of the 
spiral Eq. (10) grows quadratically with <\> [see inset of 
Fig. [3j. Then, it turns out that the asymptotic form of 
the geometric potential in terms of s is a Coulomb one^ 
[see Fig. [3]. It is then clear that the geometric potential 
corresponds to an attraction towards the point of maxi- 
mum curvature Ri n and leads to the appearance of bound 
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FIG. 4. (Color online). Behavior of the lowest eigenener- 
gies E® as a function of the quantum number n for a RUNT 
with Ri n — 3ttI and different number of rotations Nr. The 
points are the numerical results whereas the continuous line 
corresponds to the analytical behavior of the spectrum of the 
one-dimensional hydrogen Hamiltonian. The inset shows the 
behavior of the ground state binding energy as a function of 
the inner radius of the RUNT. 
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FIG. 5. (Color online). Ground state density probability as a 
function of s — Si n for a RUNT with Ri n = 3irl performing up 
to three rotations. The total length of the spiral corresponds 
to L ~ 80/ for Nr = 1,L~ 200/ for Nr = 2 and L ~ 350/ for 
Nr = 3. The inset shows the behavior of the corresponding 
eigenvalue as function of the number of rotations Nr. For 
Nr > 1 the energy of the ground state saturates at a finite 
value. 



states. The asymptotic form of the Hamiltonian Eq. ([8| 

reads 
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By restricting to the half-space s > 0, Eq. (12) is the 
Hamiltonian of a I D hy drogen atom with a 



charge" e q = h/(4\ 



quantum 
The eigenstates and the cor- 



responding eigenenergies are thus well known. However, 
in the present situation we have to meet the Dirichlet 
boundary condition at s in and s in + L. The effect of 
these boundary conditions can be captured in a two step 
process. First, the boundary condition at Si n is met by 
the infinite set of localized atomic-like states that, apart 
from a normalization constant, read 



2s 



(13) 



where U is the confluent hypergeometric function of the 
second kind and we defined the "Bohr radius" ao = 



h 2 /(r 



16/. Finally, the parameters r] n , which de- 



pend on Si n , determine the binding energies 
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Obviously, for Si n = the energy spectrum reduces to 
the usual Rydberg series (r] n = n). By increasing s^ n , 
the r]nS grow linearly with Si n meaning that the binding 
energies are inversely related to the inner radius of the 
RUNT [see inset of Fig. 

Next, we introduce the Dirichlet boundary condition 
at the outer radius of the RUNT. The atomic-like states 
do not meet this boundary condition since they do not 



vanish exactly at s out = Si n + L. However, if s out resides 
in their exponential tail, the effect of the latter boundary 
condition can be neglected. This will obviously occur for 
the lowest energy states for which s out is much larger than 
the average arclength (s). Their corresponding binding 
energies will be then accurately predicted by Eq. ( 14 ) as 
shown in Fig. |1J This is not verified for large n since the 
atomic-like states are localized over a region much larger 
than the total length of the spiral. The confinement due 
to the Dirichlet boundary conditions will dominate in the 
latter case and hence we expect the exact eigenfunctions 
to be similar to the standing waves of Eq. (11). 

By increasing the number of rotations Nr or equiv- 
alently the total length of the spiral L, one then finds 
a continuous evolution from free particle states where 
the eigenfunction is localized over the entire length L, 
to atomic-like states where the localization region is of 
the order of (s) [see Fig. |5j. Accordingly, as shown in 
the inset of Fig. |5j the eigenvalue scale with 1/L 2 in the 
free particle region saturating at the finite negative value 
given by the binding energies Eq. (14). 



Now we can determine the appearance of zero energy 
eigenstates. The asymptotic Hamiltonian Eq. (12) ad- 
mits a zero energy eigenstate which has, apart from a 
normalization constant, the following general form 



^oO) 



Ji 
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(15) 



In the equation above, J and Y indicate respectively the 
Bessel functions of the first and second kind whereas C 
is an arbitrary constant that can be fixed by requiring 
Eq. (15) to meet the Dirichlet boundary condition at the 
inner radius of the RUNT. Here it is convenient to write 
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the eigenstate Eq. (15) in terms of the azimuthal angle 
of the Archimedean spiral. Since in the (j) » 1 regime 
s ~ I (j) 2 /2 [see Eq. Hm], we find 
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(16) 



where we got rid of the Bessel functions appearing in 
Eq. ( 15 ) by taking advantage of their asymptotic expan- 
sion for large qb values 32 . From Eq. (16) it is immediately 



clear that the zero energy eigenstate meets the second 
Dirichlet boundary condition at the outer radius of the 
RUNT only for (j) out = (j)i n + 2irk with k integer and hence 
for an integer number of rotations independent of the in- 
ner radius of the RUNT. Notice that the zero energy state 
will have Nr — 1 nodes and thus will represent the Nr- 
th lowest energy state as indeed numerically found. Thus 
the number of curvature-induced bound states is equal to 
the number of windings of the Archimedean spiral. 



IV. CONCLUSIONS 

In conclusion, we have investigated theoretically the 
single particle states in a rolled-up nanotube and have 
found that the effect of the curvature results in the ap- 
pearance of atomic-like localized states. Interestingly the 
number of the bound states corresponds to the rotation 
number of the nanotube. We have also determined how 
the binding energies depend on the other relevant geo- 
metric parameters, namely, the radial super lattice con- 
stant and the typical radius of the nanotube. 
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